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Abstract 

Each platonic solid defines a single-qubit positive operator valued measure (POVM) 
by interpreting its vertices as points on the Bloch sphere. We construct simple cir- 
cuits for implementing this kind of measurements and other simple types of symmetric 
POVMs on one qubit. Each implementation consists of a discrete Fourier transform 
and some elementary quantum operations followed by an orthogonal measurement in 
the computational basis. 

1 Introduction 

A key postulate of textbook quantum mechanics is the assumption that measurements 
correspond to self-adjoint operators A in such a way that the probability of each possible 
measurement outcome or set of possible outcomes can be computed from the spectral 
projections of A. If the corresponding system Hilbert space is finite dimensional A can 
be written as A = • XjPj where Pj is the projection on the eigenspace with eigenvalue 
Xj. The probability of the outcome Xj is tr (pPj) if the system is in a state with density 
matrix p. This type of measurement is called von-Neumann measurement, orthogonal 
measurement or projector-valued measurement. 

Within the standard model of a quantum computer one can easily show that it 
is in principle possible to implement measurements for all self-adjoint operators A 
acting on the Hilbert space (C 2 )®', i.e., the state space of a quantum register with I 
qubits. Since a universal quantum computer allows the implementation of each unitary 
transformation one could perform a unitary operation U that diagonalizes A with 
respect to the computational basis and measure with respect to this basis. 

However, the description of measurements by self-adjoint operators is not gen- 
eral enough. Most general measurements are described by positive operator valued 
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measures (POVMs). A POVM is denned as follows |J. Let O be the set of possi- 
ble outcomes and S be a sigma-algebra of measurable subsets of fi. Let V be the 
set of positive operators acting on the Hilbert space 7i. Then a POVM A is a map 
A : S — > V,m i— > A m , with the following properties: 

1. For all countable families (rrij) of mutually disjoint sets rrij one has 



The probability for obtaining an outcome in the set m is given by tr(pA m ). When the 
set f2 of possible outcomes is finite or countably infinite a POVM is uniquely given 
by a family (Aj) of positive operators such that pj = tr(pAj) is the probability for 
obtaining the outcome j. We only consider POVMs with a finite set SI of outcomes. 
Furthermore, the considered POVMs have the following properties: 

1. The family (Aj) describes a single-qubit measurement, i.e., the system Hilbert 
space is C 2 . 

2. Each Aj is a rank-one operator, i.e., Aj = \^fj)(^j\. The vectors \^j) have the 
same length. They are not necessarily normalized. 

3. The operators Aj correspond to symmetric points on the Bloch sphere. The 
symmetry groups are finite subgroups of 50(3). The possible symmetry groups 
are the cyclic and dihedral groups and the symmetry groups of the platonic solids. 

These properties show that we restrict our attention to a rather specific class of 
symmetric POVMs. The symmetry is fundamental in our constructions of the circuits 
implementing the POVMs. Specifically, we choose a cyclic subgroup of the symmetry 
group corresponding to a POVM. Under the action of the cyclic group the set of points 
on the Bloch sphere decomposes into several orbits. As shown in Section |I] POVMs 
given by a single orbit can easily be implemented by a discrete Fourier transform. Since 
we have several orbits we have to use additional gates besides the Fourier transform 
to implement the POVM. This explains why the discrete Fourier transform plays a 
central role in all constructed circuits. 

The intention of this paper is to show how the symmetry of a POVM can be used 
to construct a simple circuit for implementing the POVM. To our knowledge, there 
are no considerations of the implementation of POVMs besides !2 t . The investiga- 
tion of the implementation and its complexity is motivated by the fact that there are 
examples where generalized measurements can extract more information about an un- 
known quantum state than projector- valued measurements. Symmetric POVMs may, 
for instance, be interesting when we want to distinguish between symmetric states j2j. 
Furthermore, POVMs may perform better than orthogonal measurements with respect 
to appropriate information criteria (e.g. mutual information (3) or the least square er- 
ror [1]). Here we do neither consider these "quality" criteria nor the post-measurement 
state. The post-measurement state may be relevant in order to understand information- 
disturbance trade-off relations 5 . 




j 



where the infinite sum converges in the weak operator topology. 
2. An = 1. 
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In the next section we describe the basic principles for implementing arbitrary 
POVMs. In Section we specify the correspondence of POVM operators to points on 
the Bloch sphere. Furthermore, we specify the symmetry of POVMs. In Sections0]and 
El we consider the implementation of POVMs with a cyclic or dihedral symmetry group, 
respectively. These considerations are the basis of the implementations of POVMs 
corresponding to platonic solids. The implementation of these POVMs is discussed in 
Sections EHES 

2 Orthogonal measurement of POVMs 

In this section we briefly rephrase Neumark's theorem describing the reduction of 
POVMs to orthogonal measurements .6.. This theorem allows to implement POVMs 
by performing unitary transformations on the joint system consisting of the system to 
be measured and an ancilla register. The unitary transformations are followed by an 
orthogonal measurement in the computational basis. 

Let (Aj) with j £ {l,...,n} be a POVM with corresponding Hilbert space C d 
where each Aj = \^>j)(^j\ £ <C dxd is a positive operator of rank one. Due to the 
properties of POVMs we have ^ ■ Aj = Id where Id denotes the identity matrix of 
size d. The choice of corresponding vectors Y&j) is not unique since we can multiply 
each l^j) with a phase factor that is physically irrelevant. It is therefore reasonable 
to choose the phase factors in such a way that the implementation of the POVM is 
simplified. Our constructions in Sections 4-10 implicitly make use of this. For n > d 
the vectors \^j) cannot be mutually orthogonal. As a simple example we consider a 
system with Hilbert space C 2 and the following vectors: 

i*,)=^(;), i* 2 > = U ) - d i*=> = \/K i>) • 

Here is uj := exp(— 2m/ 3) a third root of unity. We therefore have 

as POVM operators. In Section H] we consider a generalization of this POVM. 

Assuming orthogonal measurements as basic measurements, we have to extend 
the system by at least n — d dimensions to make a measurement with n different 
measurement outcomes possible. In order to simplify notation, we consider the given 
system with d dimensions as a subsystem of a system with n dimensions. Since we 
are interested in quantum circuits we have to embed the system into a qubit register. 
This can be done by assuming that the POVM consists of n = 2 l operators. Note 
that this is no loss of generality since we can extend a given POVM by an appropriate 
number of zero operators Aj = 0^ S Qdxd w j iere g d denotes the zero matrix of size 
d. This extension does not change the probability distribution of the POVM since 
Pj = tr(pOd) = for a zero operator Aj = 0d- In our example above we add the zero 
operator A^ = O2 to the three POVM operators. We obtain a POVM that can be 
implemented on a register of two qubits. 
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The basic idea of Neumark's theorem is to implement an orthogonal measurement 
(Aj) on the extended system with n dimensions that corresponds to the POVM (Aj) 
in the sense that it reproduces the correct probabilities pj. We now consider the 
construction of the orthogonal measurement {Aj). Let p G C dxd be the density matrix 
of the state to be measured. Then the state of the extended system with n dimensions 
can be written as p = p(BO n ~d S C nxn . When we write the vectors \^>j) as columns of 
the matrix 

M=(|* 1 )...|$ ri ))GC (ixn , 



the operators Aj = G C nxn are given by 

extended vectors are the columns of the matrix 



1$,) G C n . The 



M = 



l*i> 
l*i> 



G C 



that is an arbitrary unitary matrix containing M as upper part of size d x n. The 
extension of M to a unitary matrix M is always possible since the rows of M are 
orthonormal. This is guaranteed by the fact that each POVM (Aj) satisfies ^ ■ Aj = 
The probability distribution pj = tr(pAj) equals the distribution pj of the original 
POVM since 



Pj = tr [pA^j = tr (ip © n _ d ) 



= tr (Mi) =Vy 



In our example, we have 



M 



1 



or 



G C 



2x4 



A possible unitary extension M of this matrix is given by 



M 



leading to the vectors \^\) 
y/T/3(l,u 2 ,u,0) T , and |* 4 ) = (0,0,0, 1) T . 
instance, we obtain the probability 



P2 = tr 




-<4x4 



1/3(1, u,Lu 2 ,0y 
With the state p = p © 2 G 



1*3) = 
: 4x4 , for 



/ 


/ 


Pn 


Pl2 





o \ 


/ 


1 


C^ 2 


UJ 


M 


\ 




P21 


P22 












1 


UJ 2 





























1 







V 


V 











0/ 


V 











o J 


J 



tr p 



1 



P2 



for the second POVM operator. The probabilities for the other POVM operators are 
computed similarly. 

The implementation of a POVM with corresponding matrix M is obtained by the 
orthogonal measurement in the computational basis after performing the unitary trans- 
formation Aft on the system with initial state p. This unitary operation maps the vector 
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\^fj) = © \&j) to the computational basis vector \j). Therefore, the measurement 
in the computational basis after applying corresponds to the measurement in the 
basis defined by the vectors 

In summary, we are interested in constructing and implementing the matrix M' for a 
given POVM corresponding to the matrix M. In the following sections the construction 
of the matrices M' is considered for symmetric POVMs besides the decomposition 
of M' into elementary (one- and two-qubit) gates. The symmetry leads to simple 
constructions and implementations based on Fourier transforms. 



3 Symmetric POVMs on a single qubit 

As described in the previous section, the basis of the orthogonal measurement of a 
POVM with corresponding matrix M is the implementation of M'. M is a unitary 
extension of M. We can apply the algorithm in Section 4.5.1 of 7 to obtain a quantum 
circuit for M\ The algorithm decomposes the matrix M> into a product of two-level 
matrices that can be translated into a sequence of elementary gates, i.e., each gate op- 
erates on one or two qubits. In general, the constructed circuit for is of exponential 
size in the number n of POVM operators. Intuitively, some symmetry properties of 
the considered POVMs may lead to algorithms constructing smaller circuits than the 
standard algorithm that works for arbitrary unitary matrices. 

To specify the symmetry of POVMs on a single qubit using geometric concepts, we 
use the correspondence of POVMs to points on the Bloch sphere as already mentioned 
in the introduction. Usually, each point on the Bloch sphere is considered as a pure 
state. Specifically, a pure state p G C 2x2 corresponds to the point {x,y,z) T G K 3 on 
the Bloch sphere with 

x \ / tv(a x p) 
V = tr (°VP) 



where 



z J \ tr(<7 2 p) 



, 1 \ / -i \ . f I 

(Tt, = h „ , <r„ = . „ , and a z 



i o ) ' y i o r ' v o -l 



denote the Pauli spin matrices. Conversely, the point (x, y, z) T € M 3 on the Bloch 
sphere corresponds to the density matrix 

1 ( 1 + z x-iy \ ^_ ^2x2 



2 \ x + ly 1 — z 

For some special states the points on the Bloch sphere are shown in Figure ^ We 
now extend the Bloch sphere representation for states to a representation of POVM 
operators of rank one. Note that each pure state is a projection p = |\&)(^ / | of rank 
one. By rescaling a POVM operator Aj = l^,-}^^ to a density matrix we can identify 
Aj with a point on the Bloch sphere using the correspondence for pure states. 

A symmetry of a POVM can be defined by a symmetry of the corresponding points 
on the Bloch sphere. We are interested in POVMs with a finite symmetry group for the 
points on the Bloch sphere, i.e., we consider finite subgroups of SO(3) |S]. There are 



5 




two infinite families of finite subgroups, namely the cyclic groups C m and the dihedral 
groups D m for m > 2 as symmetry groups of an m-sided regular polygon (with the 
special case m = 2). Furthermore, we have the symmetry groups of the five platonic 
solids (tetrahedron, cube, octa-, dodeca-, and icosahedron). As a restriction for the 
latter symmetry groups, we assume the points on the Bloch sphere of a POVM to 
coincide with the vertices of the platonic solid corresponding to the symmetry group. 

The vertices of the regular polygons and the platonic solids depend on the orien- 
tation of the polygons and platonic solids in the Bloch sphere. Each orientation leads 
to another POVM. In order to simplify the constructions in the following sections we 
choose specific orientations of the regular polygons and platonic solids. To obtain the 
implementation of a POVM corresponding to the same polygon or platonic solid with 
another orientation, it suffices to implement a single qubit operation on the qubit to 
be measured. 



4 Cyclic groups 

The simplest finite symmetry groups of points on the Bloch sphere are the cyclic 
groups. For a fixed m > 2 we consider the rotations of an m-sided regular polygon 
with a common axis perpendicular to the face of the polygon. The rotations form a 
group that is isomorphic to the group C m = (r) with r m = 1. The implementation of 
POVMs corresponding to a single orbit of points under the action of a cyclic symmetry 
group is the basis of all constructions in the following sections. 

In principle, we can choose an arbitrary orientation of the polygon corresponding 
to the cyclic symmetry group. For simplification, we choose the face of the regular 
polygon to be perpendicular to the z-axis in the Bloch sphere. In other words, the 
z-axis is the common axis of the rotations. For instance, the 5-sided regular polygon 
is shown in Figure |5J 
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Figure 2: Points of the cyclic POVM in the xy-p\a.ne for m = 5. 



The cyclic symmetry group of the polygon is generated by the 27r/m rotation about 
the z-axis. In the Hilbert space, this rotation of the Bloch sphere corresponds to the 
matrix diag(l,cj) G C 2x2 where uj := exp(— 2iri/m) is an mth complex root of unity. 
The diagonal form of the matrix is the reason for choosing the z-axis as common 
rotation axis. If we choose the vector (1, 1) T G C 2 and consider the orbit under the 
symmetry group then we get the vectors (1, uJ) T G C 2 for j G {0, . . . , m — 1}. Other 
vectors of the Bloch sphere do not lead to POVMs or to POVMs that correspond to a 
polygon with another rotation. The latter case is discussed at the end of the previous 
section. Due to the identity 

the elements |%) = y/l/m(l, iJ~ 1 ) T G C 2 for j G {l,...,m} define a POVM with 
m = n operators on a qubit. Therefore, the unitary matrix M G <C mxm is a unitary 
extension of the matrix 

M = JI.( ) 1 ••• Li ) £C 2xm 
V m \ 1 a; ... uj m 1 J 

M corresponds to the first two rows of the discrete Fourier matrix 

771 / 1 / jk\ m ^ — rr<rnxrn 

V m V / i,fe=o 

of size m. Consequently, by considering the qubit to be measured as a subsystem of an 
m-dimensional system the implementation of the inverse Fourier transform = Fm 
leads to the probability distribution of the cyclic POVM on the qubit. 

For the construction of a quantum circuit we have to embed the system of dimension 
m into a register with I qubits. The register must have r := 2 l > m dimensions. 
Following Section [21 we extend the cyclic POVM by an appropriate number of zero 




Figure 3: The general circuit scheme for implementing the cyclic POVM (left side) and the 
circuit for implementing the cyclic POVM for m = 4 (right side). 



operators A m+ i, . . . , 4=0 2 GC 2x2 . We therefore have 

n \ 

G C 



m \ 1 to ... <jj 



TO— 1 



... \ _ „2xr 



... 



A possible unitary extension M of this matrix is given by M — F m © Ir—m 
where I r -m denotes the identity matrix of size r — m. Consequently, on a qubit register 
the cyclic POVM corresponding to the m-sided regular polygon can be implemented 
by performing the operation Aft — F m © J, — m . The circuit for implementing the cyclic 
POVM is schematically shown in Figure |3J Note that the embedding p i— > p © r _ m 
corresponds to the use of initialized ancilla qubits. 

The Fourier transform can be implemented efficiently if m is a power of two [7j. 
Furthermore, the embedding into a qubit register is straightforward since we do not 
need zero operators in this case. In summary, the cyclic POVM can be implemented 
efficiently on a qubit register if m is a power of two. For instance, the quantum circuit 
for the implementation of the cyclic POVM is shown in FigureOJfor m = 4. The circuit 
of f\ is the standard circuit for Fourier transforms |7J . Note that the first permutation 
of the qubits can be removed when we change the order of the input. 



5 Dihedral groups 

The cyclic symmetry group of an m-sided regular polygon which we considered in the 
previous section is a subgroup of the dihedral group. The dihedral group consists of all 
rotations which map the m-sided regular polygon onto itself. In contrast to the cyclic 
group we allow the rotations to have different axes. For a fixed m > 2, the dihedral 
group is isomorphic to D m = (r,s) with r m = 1, s 2 = 1, and srs = r _1 . In order to 
use the results for the cyclic groups, we consider the same orientation of the regular 
polygon as in the previous section, i.e., the face of the polygon is orthogonal to the 
z-axis. Furthermore, we assume that at least one vertex is an element of the x-axis. 
Due to this orientation, the element r corresponds to the 2ir/m rotation about the 
z-axis and the element s corresponds to the ir rotation about the x-axis. In the Hilbert 
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Figure 4: Points of the dihedral POVM with m = 5. 



space these rotations correspond to the matrices 





u 



and 



where uj := exp(— 2m /m) is a mth complex root of unity. We can define a projective 
representation of the group D m by mapping the element r £ D m to the first matrix 
and the element s £ D m to the second matrix. 

We consider the orbit of a vector under the action of the dihedral group D m . Let 
(a, j3) T G C 2 with |a| 2 + \f3\ 2 = 1. Since a global phase factor of a vector is physically 
irrelevant we assume a£l without loss of generality. Under the action of the dihedral 
group the orbit contains the vectors (a, /3u;- ? ) T and (/?, aLO^) T with j G {0, . . . , m — 1}. 
An example of the orbit is shown in Figure 0] We have at most n = 2m vectors. In the 
following, we assume that the orbit contains 2m elements. If the orbit of D m contains 
less than 2m points we have either the case that all points are on the xy-plane (and the 
POVM consists of a single orbit under the group C m ) or we have only the two points 
(1,0) T and (0, 1) T defining an orthogonal measurement. Since 

E ( ^ ) («>^') + E ( L ) ^ auJ ~ 3) = m/2 ' 



we rescale a and (3 with the factor ^Jl/m to obtain a POVM. 

We now consider the implementation of the dihedral POVM. In order to analyze 
the structure, we do not consider the embedding of the constructed system into a qubit 
register in the first place. The orbit under the action of D m breaks into two orbits under 
the action of the subgroup C m . The two orbits can be obtained by the action of C m 
on the vectors {a,[3) T and a) T . Therefore, we expect to obtain implementations of 
the dihedral POVMs which are similar to the implementations in the previous section. 
With an appropriate order of the vectors we have the matrix 
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/3 (3u 
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-<2x2m 
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For even m, this matrix can be extended to the unitary matrix 



M = Q 
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with a permutation matrix Q £ C nxn fixing the first row and mapping the (m + 2)nd 
row to the second row. For odd m, the extended matrix is similar. We only have to 
write 0cv^ m ~ 1 ^ instead of —0u^ n ~ 1 ^ in the mth row. In the last row we write /?u/ m-1 ^ 
and — ct!u/ m_1 )- 7 instead of 0co^ m ~ 1 ^ and ccu/" 1-1 ^ , respectively. In order to simplify 
notation, we mainly consider the case of even m in the following. The constructions 
for odd m are similar. 

We consider a decomposition of the matrix Q'M to obtain a decomposition of M. 
The matrix Q^M can be multiplied with I2 <8> from the right leading to 

diag(a, a, a, . . . , a, a) diag(+/3, -0, +0, . . . , +0, -0) 
diag(/3, 0,0,..., 0,0) diag(— a, +a, -a, . . . , -a, +a) 

We now embed the system with n = 2m dimensions into a qubit register. We consider 
a register with I qubits where r := 2 l > re. We replace the matrix T with the matrix T r 
of the same structure but of size r. This is done by extending each of the four diagonal 
components to a diagonal matrix in C( r / 2 ) x( ' r / 2 ) while conserving the structure. For 
instance, the matrix 
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is extended to the matrix 
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Furthermore, in the factorization T = QftMiJi (8> Fm) the matrix Q is replaced by a 
permutation matrix Q r G C rxr that fixes the first row and maps the (r/2 + 2)nd row 
to the second row. In qubit notation, this permutation matrix can be described by 
|0 . . . 0) i — > |0 . . . 0) and 1 10 . . . 01) i— > 1 00 . . . 01). This permutation can be implemented 
by an XOR-gate on the first qubit controlled by the last qubit. Other implementations 
that satisfy the two constraints are also possible. The Fourier transform F m is replaced 



by F n 



'r/2— m- 



In summary, we obtain a matrix M r that is defined by the equation 



M r = Q r T r (I 2 (8) (F m I r/2 -m)) G C 



(1) 



This matrix is a unitary extension of the matrix M corresponding to the dihedral 
POVM with some zero operators as discussed in Section |2 Our example with Tg leads 
to the matrix 
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The matrix Q$ maps the sixth row to the second row leading to the first two rows 

P P P ( 



a a a 
P Puj Puj 2 



a auj au> 



with two zero columns that do not change the POVM due to zero probability. For con- 
venience, we shift the qubits according to the mapping \x± . . . xi_\Xi) i— ► \xix\ . . . x^_i). 
We denote this permutation by R. After this reordering of qubits the matrix T r takes 
the simple form 



RT r R) = A <g) J r/4 := y/m 
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By combining this equation with Equation Q we get the factorization 
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r/4) 



Translating this equation into a quantum circuit, the decomposition of M r leads to the 
circuit scheme shown in Figure [5J The operation is decomposed as 
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Figure 5: A quantum circuit for implementing the dihedral POVM. We set a := ^fma and 
b := to simplify notation. 

corresponding to the second and third gates from the left in Figure |SJ We do not 
have to implement the permutation R explicitly if the controlled one-qubit operations 
are applied to appropriate qubit pairs. The given circuit can be slightly simplified by 
merging the first two gates from the left to a single controlled gate with the operation 

As discussed in the previous section the Fourier transform F m can be implemented 
with a polylogarithmical number of gates if m is a power of two. Consequently the 
dihedral POVM can be implemented efficiently in these cases. 

6 Tetrahedron 

The tetrahedron is the platonic solid with four faces. The symmetry group of the 
tetrahedron is isomorphic to the alternating group A4. This group consists of the 
twelve permutations of four elements with positive signum. We consider the POVM 
corresponding to the vertices of the tetrahedron in the Bloch sphere. The tetrahedron 
is shown Figure El For instance, the vertex 1 is given by the vector (-v/2/3, 0, a/1/3) t € 
M 3 . The vertices of the tetrahedron correspond to the vectors 




with a = v/(3 + V3)/6 and (3 = yj (3 - y/S)/6. The first pair of vectors corresponds 
to the vertices 1 and 2, the second pair corresponds to the vertices 3 and 4. Note the 
similarity of these vectors to the vectors 




F m ® l r/2-m 
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Figure 6: The tetrahedron with two edges perpendicular to the z-axis. 



These vectors result from the action of the dihedral group with m = 2 as considered 
in previous section with the vector (a,P) T . The factor i in the second component of 
the last two vectors of Line (J2J) results from the 7r/2 rotation about the z-axis of the 
lower edge with vertices 3 and 4 relative to the upper edge with vertices 1 and 2. This 
rotation corresponds to the matrix diag(l,i) G C 2x2 . Due to the equation 

(a, 13) + 
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we have the matrix 
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P \ g r 2x J 
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with the rescaled elements a = y (3 + \/3)/12 and /3 
can be extended to the unitary matrix 

I a a (3 (3 
a —a —(3i Pi 



(3 - v / 3)/12. This matrix 
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acting on a register of two qubits with the permutation matrix 
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The matrix Q can be implemented by an XOR-gate on the first qubit controlled by 
the second qubit. We consider the decomposition of Q^M to obtain a decomposition 
of M. After multiplying Q^M with (J 2 <g> F 2 ) G C 4x4 we have 



Q^M(I 2 ®F 2 ) = V2 
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Figure 7: The circuit for the tetrahedral POVM. We set a := \f2a and b 
notation. 



\f2(3 to simplify 



and after multiplying this matrix with diag(l, 1, 1, i) £ C 4x4 from the right we have 
the equation 



Q ] M{h 
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V2 



a (3 
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h € C 



4x4 



(3) 



The matrix diag(l, 1, 1, i) corresponds to a controlled phase gate diag(l, i) on the second 
qubit. Using Equation © we get the equation 



( 1 \ 

10 

10 

V i J 



V2 



a (3 
/3 —a 



)t 



Consequently, the circuit in Figure [7| implements the transformation for the POVM 
corresponding to the tetrahedron. 



7 Cube 

The POVM associated with a cube in the Bloch sphere is a special case of the dihedral 
POVMs considered in Section |S] with m = 4. Nevertheless, we consider the implemen- 
tation of the cubic POVM in this section since we can obtain a smaller circuit by using 
the special values of a and (3. As in Section |S] we rotate the cube in the Bloch sphere to 
obtain a face perpendicular to the z-axis. Furthermore, we can rotate the cube about 
this axis to get points corresponding to the vectors 

(/?)'( pi )'(-/?)'( -A - ) '(-«)'( ) ' ( « ) ' ( «0 g £2 

with a = \J (3 + \/3)/6 and /3 = ^/(3 — \/3)/6. The first four vectors correspond 
to vertices 1-4 in Figure El and the last four vectors correspond to vertices 5-8. For 
instance, the vertex 1 corresponds to the Bloch point (^2/3, 0, ^1/3) T G R 3 . Note 
that a and (3 are real numbers. This allows us to use a more efficient construction than 
in Section |SJ Since we have the equation 

( ; ) («. P) + -+( ) («. 00 + ( _1 ) (A -«) + ...+ ( £ ) (A -ai) = 4I 2 
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Figure 8: The cube with two faces perpendicular to the z-axis. 



the given vectors define a POVM when we rescale a and with 1/2. The matrix M 
corresponding to the POVM is given by 



M 



a a 
p pi 



a 



-(3 



a 

-pi 



P 



-a 



P 

-ai 



P 
a 



P 
ai 



g 



-.2x8 



In contrast to the construction of the dihedral POVM we use the fact that for even 
m the element —1 is in the set {l,a>, . . . ,uj m ~ 1 } where u := exp(— 2iri/m) is an m-th 
complex root of unity. Therefore, we can reorder the vectors (P, au>i) T G C 2 considered 
in SectionElto obtain a matrix M with the partial row (—a, —ato, . . . , -au m_1 ) instead 
of (a,aw, . . . ,auj m ~ 1 ). This is besides G R the second reason that allows a more 
efficient construction compared to the construction for the dihedral POVM. Using the 
equation /? = P we can extend M to the unitary matrix 



M = Q 
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acting on a register of three qubits with a permutation Q satisfying 1 000) \—* |000) and 
1 101) i — ^ |001) in qubit notation. For instance, the permutation can be implemented 
by a single XOR-gate on the first qubit controlled by the third. We now consider the 
special structure of Q^M to obtain a decomposition of M*. More precisely, the matrix 
Q^M can be written as the following tensor product 



This leads to the identity 

M f = 



a 

P 



P 



-a 



a P 
P —a 



F Q T eC 



G C 



8x8 



-8x8 
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Figure 9: The circuit implementing the cubic POVM. We set a 
notation. 



2a and b := 2/3 to simplify 




Figure 10: The octahedron with two faces perpendicular to the z-axis. 

defining the quantum circuit given in Figure EJ Compared to the general circuit in 
Section we are able to replace two controlled gates by a single uncontrolled gate. 



8 Octahedron 

The symmetry group of the octahedron is identical to the symmetry group of the cube 
since the octahedron is the dual polyhedron of the cube. The group is isomorphic to 
the symmetric group S4. This group consists of all 24 permutations of four elements. 
A simple implementation of the octahedral POVM can be obtained by the orientation 
of the octahedron as shown in Figure ^] where the upper face with vertices 1-3 and 
the lower face with vertices 4-6 are perpendicular to the z-axis. Vertex 1 corresponds 



to the real vector (y 2/3,0, 
points 1-6 are given by 



1/3) G R . The complex vectors corresponding to the 



a 

Puj 



a 



Pu 2 



P 



-a 



P 

-olu 



P 



-aw 



E C 2 



where ui := exp(— 2tti/3) is a root of unity, a = y (3 + \/3)/6 and P = y (3 — \/3)/6. 
The first three vectors correspond to the upper three vertices 1-3 of the octahedron, 
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the last three vectors to the lower three vertices 4-6. Despite the negative sign of 
the second component of the last three elements, these vectors are identical with the 
vectors 



a 



au 



au 



The latter vectors are obtained by the vector (a, 0) T under the action of the dihedral 
group D3 as discussed in Section |21 Similar to the factor i in two vectors of the 
tetrahedral POVM in Section 13 the negative sign results from the tt rotation about 
the z-axis of the lower three vertices 4-6 relative to the upper three vertices 1-3. Since 



(«,/?) + 







-au 



(0, -au) = 3I 2 



we rescale a and with the factor W 1/3 to obtain a POVM. Therefore, we have 



M 



a 



a 
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(3 (3u 0u 2 



[3 (3 
—a —au 







-au 



^2x6 



(4) 



As already discussed, we have the negative signs in Equation (J3J because the lower face 
is rotated relatively to the upper face. This is different from the cubic POVM where 
we have to reorder the operators (compared to the dihedral case) in order to get the 
negative signs in the second component of the last four vectors. To see this, we write 
these components as —a, -aw, —au 2 and —au 3 with u = i as in the dihedral case. 

We now consider the extension of M to a unitary matrix M. As in Section El we do 
not embed the system with six dimensions into a qubit register in the first place. The 
matrix M corresponds to the first two rows of the matrix 



M = Q 



where Q is a permutation matrix that fixes the first row and maps the fifth row to the 
second. Similar to the previous section, this matrix can be written as 
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M = Q[V3 



a 
—a 



ft eC 



6x6 



(5) 



We now translate the decomposition of M into a circuit. We have to embed the 
system with six dimensions into a qubit register with at least three qubits. This can be 
done by replacing the Fourier matrix F$ in Equation Q with F3 © I\ where I\ G C lxl 
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Figure 11: A circuit for implementing the octahedral POVM. We set a 
to simplify notation. 



\/3a and b := y/3/3 



denotes the identity matrix of size one. This replacement leads to the matrix 
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where Q% E C 8x8 is a permutation matrix that fixes the first row and maps the sixth 
row to the second row. In qubit notation, these constraints are given by |000) *— > |000) 
and 1 101) i-^ 1 001). For instance, this transformation can be implemented by an XOR- 
gate on the first qubit controlled by the last qubit. If we restrict Ms to the first two 
rows we get the matrix 



a a a 

010 0LO 2 







-a: 





-au 







-aw 



corresponding to the desired POVM. The POVM operator corresponding to the fourth 
and eighth column is O2 £ C 2x2 leading to a zero probability for all states p E C 2x2 . 
In summary, we have the equation 



.W. 







-a 



fJ e h 



qI 



for the implementation of the octahedral POVM. This equation corresponds to the 
circuit shown in Figure ITT1 



9 Dodecahedron 

The dodecahedron is the platonic solid with twelve faces and twenty vertices. The 
symmetry group of the dodecahedron is isomorphic to the alternating group A5. This 
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Figure 12: The dodecahedron with two faces perpendicular to the z-axis. 

group contains the sixty permutations of five elements with positive signum. The 
dodecahedron is shown in Figure IT21 The upper face with vertices 1-5 and the lower 
face with vertices 6-10 are perpendicular to the z-axis. The point 

corresponds to vertex 1. This orientation of the dodecahedron in the Bloch sphere 
leads to a simple construction of the dodecahedral POVM. With to := exp(— 2iri/5), 
the points on the Bloch sphere correspond to the complex vectors 

( £ ) ■ ( -L )•{&)•{ -v ) £ c2 (6) 

where j G {0, . . . ,4}. The vectors (a,(3uj : ') T correspond to the points 1-5, the vectors 
(/?, — au}i) T to 6-10, the vectors { r y,5uji) T to 11-15, and the vectors (5,^uji) T to the 
points 16-20. The parameters a, /?, 7, and 5 are defined as follows: 

« = ^ + ^/75 + 30v5, /5=^-^/75 + 30a/5 

and 

Due to the equation 

(^)(«,/5) + ...+ ( 7 ^ 4 )(5,7O = 10/ 2 

we rescale the elements a, /3, 7, 5 with the factor yT/10 to obtain a POVM. In contrast 
to the constructions of Sections EHH1 the points on the Bloch sphere decompose into 
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four different orbits under the rotation about the z-axis. Note that there are two 
pairs of orbits. In Line (jHJ) the vectors (a, j3u^) T and (P,—au}i) T are similar to the 
vectors (a, /3a^) T and (P,aui : ') T . The latter vectors are the orbit of (a,P) T under the 
dihedral group with m = 5 as considered in Section As in the previous section a it 
rotation of one orbit relative to the other orbit causes the negative sign of the elements 



(/V 



-OLD- 



Analogously, the second orbit under is defined by the third and fourth 



type of vectors in Line ©• In summary, the vertices of the dodecahedron correspond 
to two orbits under the dihedral group -D5 with a ir rotation about the z-axis of some 
points on each orbit. Consequently, we can expect to use a similar construction as in 
the previous sections. 

We now consider the construction of the circuit for implementing the dodecahedral 
POVM. For convenience, we do not embed the system into a qubit register in the first 
place. We have the matrix 
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This matrix corresponds to the first and second row of the unitary matrix M defined 
by the equation 

M = Q {A ® F 5 ) E C 20x20 , (7) 

where Q £ c 20x2 o j s a permutation matrix that fixes the first row and maps the seventh 
row to the second row. The matrix A is defined by 



A = V5 



Now, we want to embed the extended system into a register with five qubits. Similar 
to the construction for the octahedron in Section |SJ we can do this by replacing the 
matrix F§ in Equation (JJJ) by the matrix {F^@Iz) £ C 8xS where ^3 denotes the identity 
matrix of size three. The matrix Q is replaced by a permutation matrix Q32 E C 32x32 
that satisfies 1 00000) i-> 1 00000) and 1 01001) i-> 1 00001) in the qubit notation. This 
permutation can be implemented by an XOR-operation on the second qubit controlled 
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by the last qubit. In summary, the matrix M^ 2 is defined by the equation 
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The corresponding circuit is shown in Figure 1131 Note that the matrix A' can be 
written as product A* = (h © (— <t z )) (h ® B) R (I2 <S> C) with the matrices 
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Figure 13: A circuit for implementing the dodecahedral POVM. In order to simplify notation, 
the elements ±1 of the gates on the first qubit represent ± a/1/2. 



The matrix R is the product 
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In Figure 1131 the latter two matrices correspond to the two operations on the first 
qubit which are controlled by the second qubit. 



10 Icosahedron 

The icosahedron is the dual polyhedron of the dodecahedron. Consequently, the sym- 
metry groups of both platonic solids are identical. We assume the specific orientation 
of the icosahedron as shown in Figure ITU to obtain a simple construction of the icosa- 
hedral POVM. The upper face with vertices 1-3 and the lower face with vertices 4-6 
are perpendicular to the z-axis. Vertex 1 is given by the vector 



10 - 2y/5 



15 



0. 



5 + 2^/5 



15 



The vertices of the icosahedron in the Bloch sphere correspond to the complex vectors 



a 



-auj J 



7 _ 



(8) 
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Figure 14: The icosahedron with two faces perpendicular to the z-axis. 



with j £ {0, 1, 2}, ijj := exp(— 2ni/3) and 
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1 + ^/75 + 30^5 
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75 + 30\/5 



and 
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75 - 30\/5 , 5 



2 30 

The vectors in Line © with j = correspond to the vertices 1, 4, 7 and 10 in the 
given order. As in the case of the dodecahedron we have four orbits under the rotations 
about the z-axis. Therefore, we can expect that a similar construction as in the previous 
section is possible. Due to the identity 



we have the matrix 
M = 



a a 
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where we rescale a, 0, 7 and 5 with the factor yT/6. The matrix M consists of the 
first and second row of the matrix 



M = Q (A F 3 ) € 



U2xl2 



(9) 



where Q £ c 12x12 j s a permutation matrix that fixes the first row and maps the fifth 
row to the second. Similar to the previous section, the matrix A is given by 



A = V3 
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Figure 15: A circuit for implementing the icosahedral POVM. In order to simplify notation, 
the elements ±1 of the gates on the first qubit represent ± a/1/2. 



The embedding into a register with four qubits works analogously to the previous 
section. We replace F3 in Equation © by (F3 © I\) £ C 4x4 where I\ denotes the 
identity matrix of size one. The matrix Q is replaced by the matrix Q\§ £ (£ 16x16 that 
can be described as |0000) i-> |0000) and |0101) h-> |0001) in the qubit notation. This 
permutation can be implemented by an XOR-operation on the second qubit controlled 
by the last qubit. Therefore, Equation (jJJJ) translates into 
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The circuit corresponding to this decomposition of MJ 6 is given in Figure 1151 The 
matrix can be translated into single- and two-qubit gates as shown in the previous 
section. In this translation we have to replace the constants u± and v± with 



u± 
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To 
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2± 



50 ± 5^10(5 + y/E) and v± 

11 Conclusions 



We have shown that all POVMs given by the vertices of platonic solids can be imple- 
mented using a discrete Fourier transform and a few other operations. The algorithms 
use the symmetry of the POVMs. A common feature of all constructions is the par- 
tition of the POVM operators into orbits under the action of a cyclic group. Since 
the Fourier transform allows to implement POVMs associated with an orbit under a 
cyclic group it is an essential part of all circuits. For most POVMs corresponding to a 
platonic solid, a tensor product of a Fourier transform and a specific low-dimensional 
matrix is a central building block of the circuit. The low-dimensional matrix represents 
in some sense the relations between the orbits. 

The implementation of non-symmetric POVMs seems to be a non-trivial task. It 
would, for instance, be interesting to know which POVMs can be implemented effi- 
ciently, i.e., with a number of elementary gates which grows only polynomially in the 
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number n of POVM-operators. For the symmetric POVMs considered in this paper 
the question of efficiency makes only sense for the cyclic and dihedral POVMs since the 
size of the other POVMs is fixed. For n = 2 l the complexity of the Fourier transform 
F n is only polynomial in I. Therefore, the complexity of the circuits for the cyclic and 
dihedral POVMs grows only polylogarithmically in n. 

The question of the efficiency of read-out mechanisms for a single bit has no counter- 
part in classical computer science. Complexity issues in quantum information theory 
deal not necessarily with the complexity of computational problems. They are also 
interesting in the context of measurements or state preparation procedures. However, 
there are some connections between a complexity theory of these non-computational 
quantum control problems and computational problems [HI El- Connections between 
the complexity of POVM measurements and other complexity issues may be subject 
of further research. 

The authors acknowledge helpful discussions with M. Grassl and M. Rotteler. M. 
Rotteler brought the problem of implementing symmetric POVMs to our attention. 
This work was supported by grants of the BMBF project 01/BB01B. 
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